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To extract information on hadron production dynamics in the ultrarelativistic heavy ion 
collision, the space-time structure of the hadron source has been measured using Hanbury 
Brown and Twiss interferometry. We study the distortion of the source images due to the 
effect of a final state interaction. We describe the interaction, taking place during penetrating 
through a cloud formed by evaporating particles, in terms of a one-body mean field potential 
localized in the vicinity of the source region. By adopting the semiclassical method, the 
modification of the propagation of an emitted particle is examined. In analogy to the optical 
model applied to nuclear reactions, our phenomenological model has an imaginary part of the 
potential, which describes the absorption in the cloud. In this work, we focus on the pion 
interferometry and mean field interaction obtained using a phenomenological mr forward 
scattering amplitude in the elastic channels. The p-wave scattering with p meson resonance 
leads to an attractive mean field interaction, and the presence of the absorptive part is 
mainly attributed to the formation of this resonance. We also incorporate a simple time 
dependence of the potential reflecting the dynamics of the evaporating source. Using the 
obtained potential, we examine how and to what extent the so-called HBT Gaussian radius 
is varied by the modification of the propagation. 

§1. Introduction 

The non-Abelian nature of quantum chrome-dynamics manifests in the asymp- 
totic freedom of quarks and gluons at high energy or short distance. While the 
colored quanta are, in ordinary circumstances, confined in color-neutral hadrons, 
they are liberated in extreme environments at high energy density. The study by 
the ultrarelativistic heavy ion collision enables us to explore the theoretically pre- 
dicted new state of matter called the quark-gluon plasma. 

The hydro dynamical simulations at the RHIC energy successfully reproduce the 
single-particle spectra by taking into account the radial flow in the central collisions 
and the elliptic flow in the peripheral collisions, which is regarded as an indication of 
the formation of the locally thermalized state. To extract the space-time structure of 
the matter, the two-particle momentum distribution of mesons has been measured to 
apply the technique called Hanbury Brown and Twiss (HBT) interferometry. This 
method was originally developed in microwave astronomy to measure the stellar 
radiiP The interference of a photon pair found in the observation was stimulating 
enough to play a seminal role in developing quantum optics. 2 ' A similar interference 
of an identical particle pair was found in the product of pp collisions, 3 ' and has been 
utilized to measure the space-time extension of the hadron source created in the 
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collider experiments P 1 '^ 

In the ultrarelativistic heavy ion collision, the space-time structure has been 
measured exclusively using the two-particle spectrum of identical pions. To extract 
the dynamics of the source, an elaborate choice of the coordinate system has been 
employed pl>l£i which refers to the motion of a pair. One of the axes is oriented in the 
direction of the incident beams. The other two in the transverse plane are spanned 
by the direction of the averaged transverse momentum of the pair and the direction 
perpendicular to it. They are called the longitudinal, outward and sideward direc- 
tions, respectively. An advantage of using this coordinate system is that the outward 
radius reflects the temporal structure of the source as well as the spatial one, owing 
to the finite interval of the emission time of a pair. The sideward radius is separated 
from the temporal structure, and measures a bare spatial extension. The observation 
of a pronounced difference in the two radii on the transverse plane is expected to 
provide a hadronic signature for the realization of the new state of matter.^"® Once 
the quark-gluon plasma is formed, one expects that the time evolution of the system 
becomes slower in the phase transition (or possibly, a crossover) regime, which re- 
flects in the long lifetime of the hadron source, and possibly the prolonged emission 
time. 

In prior to the RHIC experiments, the dynamics of the collective expansion was 
simulated with hydro dynamical models assuming the formation of the quark-gluon 
plasmaP^ The results indicate a prolonged lifetime of the matter, and a considerable 
difference in the radii, i? out /i? s id e — 2. More elaborate models also showed the large 
ratios, R OVL t/Rmdc > 1, respectively.^ - ^ While these models precisely reproduce 
many classes of observables, especially the elliptic flow, the measured radii unexpect- 
edly show a systematic deviation from the theoretical predictions EUSJ which has 
been called the "RHIC HBT puzzle". This work is motivated by this longstanding 
issue, and we investigate a possible origin of the deviation q 

Recently, we have examined the effects of a final state interaction, which has 
not been taken into account in the hydro dynamical simulations. 171 It has been 
commonly assumed in the hydrodynamical modeling that the frequent interactions 
among particles, which maintain thermal equilibrium of the system, become ineffec- 
tive abruptly on a certain hypersurface, and that the emitted particles propagate 
freely toward the detector after that moment. However, a more realistic picture is 
that the multiple interactions take place in the vicinity of the hypersurface while an 
escaping particle traverses the cloud formed by other evaporating particles. In an 
attempt to describe these interactions in the freeze-out stage, we have introduced 
a mean field interaction into the framework of HBT inter ferometryj^'^ assuming 
that the other conventional approximations are still valid. Our model contains an 
imaginary part of the mean field potential in analogy to the optical model used in 
nuclear reactions, which describes the absorption in the cloud. Focusing on the side- 
ward radius underestimated by the hydrodynamical simulations, we have examined, 
within a schematic model of static potential, how the modification of the propagation 



*' On the experimental side, a technique called the imaging method has been developed to avoid 
assuming the Gaussian parametrization of the source profile P^ 1 ^ 
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after the emission reflects in the HBT images extracted from the asymptotic form 
of the two-body momentum distribution. We found that the images are distorted 
depending on the sign of the real part of the mean field potential: an attractive 
interaction stretches the sideward extension, while a repulsive one shrinks it. The 
absorptive part of the mean field potential acts to modify the images in a similar way 
to the effect of attractive interaction by cutting off the particles emitted in the deep 
interior and backside region of the source. A detailed account of these results is given 
in £[2] and We note that Cramer et al. and Pratt also examined independently 
but using the different frameworks the possible distortion of the images due to the 
mean field interactionPSJ'^J 

These results, obtained in the previous work, imply an improvement of the HBT 
puzzle by taking into account the cooperative effects of the attractive interaction 
and the absorption. They motivate us to pursue this possibility more seriously. A 
suitable framework to simulate the effects of the mean field interaction may be the 
microscopic transport approach, in which the effects are incorporated in the Vlasov 
termPD It will provide a consistent description of the dynamics in the hadron phase 
with the mean field interaction until the last collision, and describe its effects caused 
throughout the dynamics. Owing to the microscopic approach,^ we obtain the 
space-time distribution of the last collision point without parameter fitting for the 
sake of the freeze-out criterion, to which HBT interferometry is possibly sensitive. To 
the best of our knowledge, the effect of the Vlasov term in the numerical simulation 
has not been examined. 

To illustrate the effects of the mean field interaction more schematically, we 
however examine its effects in the very last stage of the freeze-out process, assuming 
the hydrodynamical picture. We estimate the magnitude of the effects focusing on 
the mean field interaction after the sudden thermal "freeze-out" Although the 
system is rather dilute after the "freeze-out" , the effect of the mean field interaction 
possibly remains in the vicinity of the freeze-out hypersurface, since the extension of 
the freeze-out hypersurface obtained from the hydrodynamical simulation is much 
smaller than that of the distribution of the last collision point observed in the mi- 
croscopic approach.^ Our attempt would imply more efficient effects in the denser 
system before the "freeze-out" . Our results suggest that the effect of the mean field 
interaction should be incorporated and examined in consistent description of the dy- 
namics, e.g., the transport approach mentioned above. In our current picture, HBT 
interferometry extracts the distorted profile of the freeze-out hypersurface due to the 
modification of the propagation after the emission. 

In §4.11 we construct a phenomenological mean field potential to improve the 
previous work and to embark on an application to the realistic dynamics in the 
ultrarelativistic heavy ion collision. To study the modification of the Gaussian ra- 
dius parameters, we should address the time dependence of the mean field potential 

*' The "freeze-out" indicates the sudden freeze-out involved in the hydrodynamical modeling, 
and we study the effects of the mean field interaction after the "freeze-out" which is neglected in the 
conventional models: therefore, it causes the momentum shift of the spectrum after the "freeze-out" . 
In the following, the terminologies, freeze-out hypersurface and freeze-out time, are used in the same 
sense. 
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and clarify whether it is attractive or repulsive. Because the mean field potential 
diminishes in the expansion of the matter, we examine a time-dependent potential 
within simple space-time geometry of the evaporating source, which enables us to 
avoid overestimating its effect on the observables. The sign of the real part of the 
potential is determined from a microscopic view, assuming that the two-body scat- 
tering of the pion, which has the most abundant population, dominates the origin of 
the mean field interaction. We utilize the phenomenological analysis of the elastic 
forward tttt scattering amplitudeP^ The magnitude of the mean field interaction 
and the sign of the real part are not obvious so far. They have not been studied 
in detail incorporating the momentum dependence and isospin channel dependence 
of the elementary scattering. A relevant feature of the low-energy tttt scattering in 
hundreds of MeV is the formation of p meson resonance, which reflects in both the 
real and imaginary parts of the mean field potential. Owing to the attractive regime 
below the mass of p meson resonance, we actually obtain an attractive mean field 
interaction for the pion carrying momentum in some hundreds of MeV. The effect 
of this attractive interaction is in favor of an improvement of the HBT puzzle, as 
shown in preceding works .HDfiUJ^ID 

Using the obtained phenomenological potential, we examine the modification of 
the Gaussian radius parameters in < j4.21 We find that the deviation of the sideward 
radius and its ratio to the outward radius are improved by the attractive mean field 
interaction in the low-momentum regime. The magnitude of the effect is, however, 
not strong enough to resolve the problem by the sole effect of the mean field in- 
teraction after the "freeze-out" , since the density of the pion in the exterior of the 
freeze-out hypersurface is not sufficiently large, owing to the dynamics of the evap- 
orating source. Nevertheless, the magnitude of the effect is on the same order of 
the individual effect obtained by some upgrades of the hydrodynamical modeling.^ 
Our result suggests that the consistent description of the dynamics should be exam- 
ined incorporating the mean field interaction. Section [5] is devoted to the concluding 
remarks with some brief comments on the related works. 

§2. Distortion of the HBT images by the mean field interaction 

In this section, we show how the modification of the propagation by the one-body 
mean field potential is incorporated in the framework of HBT interferometry. Along 
the detailed account of the preceding workf^ we extend it in the two important 
points. First, we consider the time evolution of the source, whereas the previous 
work is limited to the static case. This extension is required to examine the time- 
dependent mean field potential, and the effect of the temporal structure reflected in 
the difference between the outward and sideward extensions. Second, we incorporate 
the relativity in evaluating the modification of the propagation, since the pion created 
in the ultrarelativistic heavy ion collision carries high momentum in the relativistic 
regime. 
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2.1. The correlation function by the density matrix formalism 

The correlation function C in HBT interferometry is denned by 

where k\ and k 2 are the momenta of a particle pair. Pi(k) is the probability of 
detecting a particle and P 2 {k\ ,k 2 ) is the probability of the simultaneous detection 
of a pair. If we observe an independent pair, the joint probability is factorized 
into the form, P 2 (ki,k 2 ) = Pi(k\)Pi(k 2 ), and the correlation function is given by 
C = 1, independent of momenta. The effect of the interference of a pair is found 
in a deviation of the correlation function from unity. Even in the absence of the 
interactions, we find a correlation between an identical boson (fermion) pair due to 
the Bose-Einstein (Fermi-Dirac) statistics P The correlation function of an identical 
boson pair deviates in the upward from unity owing to the symmetrization with 
respect to the labels of particles, as we show in the following. In contrast, for an 
identical fermion pair, the correlation function deviates in the downward due to 
the minus sign appearing in the antisymmetrization. HBT interferometry measures 
the space-time structure of the source utilizing such a correlation originated in the 
statistics in quantum mechanics. However, the abundant particle production in the 
ultrarelativistic heavy ion collision causes the final state interactions inherent in the 
particle species. They distort the pure statistical correlations and the extracted 
HBT images. In this work, we focus on the correlation of pion pair and a final state 
interaction described on the basis of a mean field interaction. A one-body potential 
describes the interaction between a pion and the rest of the evaporating pions. 

First, we recall how the correlation function C is related to the source distribu- 
tion function in the free case. We begin with the definition of the probabilities given 
by the statistical averages at the large time t a — > 00, 

Pi(fe) := lim (a[{t a )a k (t a )) , (2-2) 

ta-*00 

P 2 (k 1 ,k 2 ) := lim {al{t a )a\{t a )a k2 {t a )a kl {t a )) , (2-3) 

where the equal time commutation relation is imposed on the creation and annihi- 
lation operator, 

a p (t),al,(t)]=S(p- P ') . (2-4) 

The statistical average of an operator O is described using the density matrix g, 
taking the trace with any complete set of bases: 

(0(t))=Tr[gO(t)} ■ (2-5) 

In the definitions of the probabilities in Eqs. (|2-2p and (|2-3|) . we take the trace using 
the complete multiparticle state in momentum space, 

1 = X) ~ [dPl - fdPnlPl-- Pn){Pl ■■■Pn\ (2-6) 
n=0 n ' ^ ^ 
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where the first term (n = 0) is the vacuum term, |0)(0|. The probabilities are then 
given by the matrix element of the reduced one-particle density matrix p(t) and that 
of the reduced two-particle density matrix fait), respectively. We obtain them by 
taking the partial traces of the density matrix g, 

Pi(fc) = lim V"-: [dpi- - fdp n (px---p n k\g(ta)\pi---p n k) 

:= lim (k\ P (t a )\k) (2-7) 

t a — >OC 

P 2 (k u k 2 ) = lim y~]—: [dpi--- [dp n (pi---p n k 1 k2\g(t a )\p 1 ---p n k 1 k 2 ) 

:= lim (fcifca|^(t a )|fcifca> (2-8) 

t a — >oo 

where we have switched to the Schrddinger picture. Owing to the commutation 
relation (|2-4p and the definition of P 2 (ki, &2)> we find the symmetry for the exchange 
of momenta, 

{kik 2 \p2{ta)\kik 2 ) = (k2ki\p 2 {t a )\k 1 k 2 } = (kik 2 \p 2 (t a )\k 2 ki) , (2-9) 

and also for the two-time exchanges, 

(kik 2 \p 2 (t a )\kik 2 ) = (k 2 ki\p 2 (t a )\k 2 ki) . (2-10) 

In the asymptotic region where the interactions are absent, the two-particle density 
matrix element in Eq. (j2-8|) is factorized into the products of the one-particle density 
matrix elements in Eq. (|2-7|) . but maintaining the relations in Eqs. (|2-9p and (|2-10|) . 

as 

P 2 (k h k 2 ) = lim [(fei|p(t )|fei)(fe 2 |p(t )|fc 2 ) + (fei|p(t a )|fe 2 )(fe 2 | / 5(t )|fei)] 

= P 1 (k 1 )P 1 (k 2 ) + lim |(fei|/5(i a )|fe 2 )| 2 . (2-11) 

In this form, we find the correlation of the pair represented by the off-diagonal matrix 
element. The interference term is semipositive definite. 

We briefly comment on the case of an identical fermion pair. Instead of the 
commutation relation (|2-4|) . we impose the anticommutation relation to the creation 
and the annihilation operator, 

{a p (t),4,(f)} =8[p-p') . (2-12) 

We then find the relations similar to Eq. (|2-9p . but with opposite signs: 

(kik 2 \p 2 (t a )\kik 2 ) = -{k 2 ki\p 2 (t a )\kik 2 ) = - {kik 2 \p 2 (t a )\k 2 ki) . (2-13) 

On the other hand, exchanging the momenta twice, we still find the same relation 
as in Eq. (I2-I0p without the sign flip, 

{kik 2 \p 2 {t a )\kik 2 ) = (k 2 ki\p 2 (t a )\k 2 ki) . (2-14) 
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If we factorize the two-particle density matrix element in Eq. (|2-8p with the relations 
in Eqs. (|2- 13|) and (|2- 14|) being preserved, we obtain the interference term for an 
identical fermion pair, 

P 2 (k 1 ,k 2 ) = lim [(ki\p(t a )\k 1 ){k 2 \p(Q\k2) ~ (ki\p(t a )\k 2 }{k 2 \p(t a )\k 1 )] 

t a ->00 

= P 1 (k 1 )P 1 (k 2 ) - lim |(fei|,S(i a )|fc 2 )| 2 • (2-15) 

The negative sign in the interference term reflects the statistics of the fermion pair. 

Considering pion pair, the objects we have to examine are the matrix elements 
of the one-particle reduced matrix in Eqs. (|2-7p and (|2-11|) . Keeping the Schrodinger 
picture, the effects of the interaction after the decoupling from the source are included 
in the time evolution of the reduced one-particle density matrix. By using the time 
evolution operator U(t a ,to), the density matrix is represented by 

p{ta) = U{ta, to)p(t )UHt a ,t ) , (2-16) 

where to is the decoupling time when a particle is emitted from the source, and is 
regarded as the initial time of the propagation toward the asymptotic region. If we 
assume the free propagation after the decoupling, Eq. ()2-16p is explicitly written as 

p{t a ) = e -*Ho(ta-to) p ( to yHo(t a -t ) ^ ( 247 ) 

where Hq is the one-particle free Hamiltonian. 

Inserting the relation (|2T6p and the complete coordinate bases into Eqs. (|2-7p 
and (|2-11|) . we have 

Pi(k) = lim /da;i/da;2(fe|^(ta,to)ki)(«i|p(to)|a:2)(aJ2|^ t (ta,to)|fe) , (2-18) 

ta-K» J J 

P 2 (k u k 2 ) = P 1 {k 1 )P 1 (k 2 ) + lim fdx 1 fdx 2 {k l \U{t a M\xi) 

t a ^OD J J 

x ( Xl \p(t )\x 2 )(x 2 \U^t a ,to)\k 2 )\ 2 . (2-19) 

The above expressions are described using two pieces: the density matrix element at 
the decoupling time, p(to, x\, x 2 ), and the amplitude of the propagation initiating 
at x and terminating with the asymptotic momentum k, 

<p k (t , x- t a ) = (k\U(t a , t )\x) . (2-20) 

By using a set of coordinates defined by x = (x\ + x 2 )/2 and r = x\ — x 2 , the 
Wigner function f(to] X ,p) is given by 

f(t ,x,p) = Jdre- ip - r p(t ;x + ^,x-~) , (2-21) 
which is the quantum analog of the phase space distribution function at t = to . 
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In advance to embarking on the interacting case, we derive the familiar form of 
the probabilities under the absence of final state interactions. Without the interac- 
tion, the amplitude in Eq. (|2-2U|) is simply given by a plane wave 



(p k (t ,x;t a ) = e 



-iE k (t a - 



t ) e ~ik-x 



where the energy E k is the eigenvalue of the free Hamiltonian, Ho\k) = E k \k). 
Substituting the above plane wave and the Wigner function defined in Eq. (|2-2ip . 
we obtain the expression of P\(k), 



P\{k) = jdxjdp 6(p- k)f(t ,x,p) = Jdxf(t ,x,k) , 



(2-22) 



where the delta function, in the intermediate of the above equation, represents the 
momentum conservation in the free motion. For the static source or the instan- 
taneous emission at t = to, the momentum distribution at to is preserved in the 
asymptotic region, if the interaction is absent during the propagation. The decou- 
pling time in general depends on the location of the emission point, to(x), e.g., for 
our interest, owing to the expansion of the source created in the ultrarelativistic 
heavy ion collision. In this case, Pi(k) is given by collecting the contribution of 
the particles emitted on the hypersurface in the space-time, and we then obtain an 
extended form 



Pi(k) = Jd 4 x 5(t-t (x)) f(t,x,k) 
d x S(x, k) , 



(2-23) 



where, motivated by the hydrodynamical model, we have defined the source function 
on the freeze-out hypersurface, S(x, k) := 5(t — to(x))f(t, x, k). 

Performing a similar procedure for P 2 (ki, ^2) m Eq. (|2-19p . we find a simple 
form of the interference term given by the Fourier transform of the source function, 



P 2 (k 1 ,k 2 ) = P 1 (k 1 )P 1 (k 2 ) 



d 4 x S(x, k) e 



I <jX 



(2-24) 



where the momenta appearing in the above relation are the averaged momentum of 
the pair, = (k^ + k^jl, and the relative momentum, = k 1 ^ — k%- Note that 
we adopt a convention of the metric such that the inner product of four vectors is 
given by AB = A°B° — A ■ B. We then obtain the familiar form of the correlation 
function, 



C(k,q) = l + V 2 (k,q) 



d 4 x S(x, k) e 



iqx 



(2-25) 



where the normalization i](k, q) is given by 



rT 2 (fc,<7) = Pi(fci)Pi(fc 2 ) ^ P?(k) • 



(2-26) 
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The approximation in the last equality is valid for the small q regime in which 
the correlation arises. In the limit of vanishing relative momentum, the correlation 
function behaves as C(fc,0) = 2 for chaotic source, and it has the asymptotic form 
C(fc,oo) = 1. 

Note that the relative momentum has only three independent components, 
and a dependent one. Owing to the on-shell conditions of the asymptotic momenta, 



E k 2 



k\ + m 2 
k 2 + m 2 



we have a relation among the components, 

n k 



— ■ q = v k ■ q , 



(2-27) 



where the velocity v k is defined by the averaged momentum, Vk = k/Ek- The 
zeroth component of the averaged momentum is approximated to satisfy the on- 
shell condition, and is denoted by k° ~ E k = \/k 2 + m 2 . Owing to the above 
constraint, what we obtain by measuring the correlation function in Eq. (|2-25p is 
thus not the four-dimensional image of the source, but a sort of projection of the 
image on the three dimensions. Using Eq. (|2-27p . we conventionally eliminate the 
zeroth component of the relative momentum. The correlation function in Eq. (|2-25|) 
is then rewritten as 



d 4 x S(x, k) e' 



-iq(x-v h t) 



C(k,q) = l + V 2 (k,q) 
= l + 7 1 2 (k,q) 



where, on the second line, we have defined the shifted coordinate, x' 
this choice, we observe the image of a "static" source, 



dx' dt S(t,x' + v k t,k) e- iq - a 



S(x',k) = Jdt S(t,x' + v k t,k) , 



(2-28) 

(2-29) 
x — v^t. In 

(2-30) 



where the explicit time dependence of the original source function and the implicit 
dependence in the shift of the spatial coordinate, v k t, are integrated out. The 
emission at t on x is embedded in the three-dimensional image as if it occurs on 
the shifted position, x' = x — v k t. As long as we measure the three-dimensional 
image by HBT interferometry, these two emissions cannot be distinguished.^ As we 
mentioned in the Introduction, an elaborate choice of the Cartesian coordinate, with 
the outward axis being parallel to k, enables us to separate the effect of the finite 
duration of the emission. In this frame, the velocity is represented in components as 
Vk = 0,ui ong ), and the shift of the coordinate is given by 



''out " 
C side : 
C long 



Xout ~V±t, 
' -^side > 

- 2- long ^long^ i 



(2-31) 
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where the subscripts stand for the outward, sideward and longitudinal directions, 
respectively. Roughly, the outward is parallel to the line of sight, and the temporal 
interval of the emissions reflects in the apparent depth of the measured image. 

In the conventional analyses, the spatial Gaussian profile of the source has been 
assumed. For a static source, the source function in Eq. (|2-30p is assumed to be 

S(x,k) = F T exp {-^- - - , (2-32) 

where Vr = Jdt is the temporal volume, which is canceled in the correlation function 
owing to the same factors in P2(&i ; &2) and Pi(fci)Pi(fc2)- To obtain a set of useful 
formulae for the variances, we consider the second moments of the coordinate x, 
regarding the correlation function in Eq. (|2-28p as a generating function. Comparing 
the cases with and without the Gaussian parametrization, we obtain 

v±t) 2 ) , 

(2-33) 



R 2 — 
-"-out — 


((^out - 


P2 _ 

-"-side — 


inside) ' 


P 2 _ 
-"-long ~ 


: (^long) 



where, as in the conventional analysis, we have chosen the frame in which we have 
fiong = as well as t> s ide = 0. The expectation value (&) and the deviation x from 
the center of the distribution are defined by 

x = x-(x) , (2-34) 
f G S(x, k) d 4 x 



JS(x,k) 



(2-35) 



In the last part of the next section, we find a modification of the formulae in Eq. 
(|2-33p to Eq. ()2-35p by the mean field interaction. 

2.2. The effects of mean field interaction 

In the previous section, we have recalled how the Wigner function, or the density 
matrix, is related to the correlation function under the conventional approximations: 
the factorization of the two-particle reduced density matrix and the absence of the 
final state interactions. Note that the plane waves of the free propagations result 
in the kernel of the Fourier transform in Eq. ()2-25p . In this section, we evaluate 
the one-body amplitude, ipk(to,x;t a ) = {k\e~ lH<yta ~ to ^\x) , of which time evolution 
is given by the one-body Hamiltonian with the mean field interaction, and examine 
how the HBT image extracted by the integral transform is distorted. We consider a 
nonrelativistic Hamiltonian, 

k = + V{X) ' (2 ' 36) 

with a central scalar potential V{x) that is assumed to be generated in the central 
collisions. While the nonrelativistic Hamiltonian is sufficient for heavy particles, 
HBT interferometry in heavy ion collision often uses the spectrum of the pion in 
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some hundreds MeV, or even higher momentum regime. In this regime, the am- 
plitude of the lightest hadron with m ~ 140 MeV should be evaluated with the 
relativistic theory. To improve this point, we extend the formalism using the proper 
time formalismP^ 1 This extension in the evaluation of the amplitude is carried out 
in parallel to the nonrelativistic case, and the detailed account of this point is given 
in Appendix A. Here, we consider the nonrelativistic Hamiltonian in Eq. (|2-36)) for 
simplicity. 

To examine the modification of the propagation, we adopt a semiclassical ap- 
proach, applying the stationary phase approximation to the path integral form of 
the amplitude. 17 ' The standard prescription of the path integral is given in the form 
that represents the transition amplitude of a particle from an initial position x to a 
final position x a . By evaluating it using the saddle point method, the action S c of 
the classical trajectory provides the amplitude, 

(x a \e- iA ^-^ \x) ~ A(x a , t a - x, t ) jMxa,t a ;*,t ) > (2 . 37) 

where the classical action and trajectory are specified by the boundary conditions, 
x(to) and x a {t a ). Recall that (fk(to, x; t a ) is the amplitude to detect a particle with 
the asymptotic momentum k at a sufficiently large time, t a . This amplitude is, 
however, proportional to the amplitude given by Eq. (|2-37p as long as the classical 
trajectory is uniquely specified by the boundary conditions x at to and k at t a — > oo. 
We prove this point in the following. 

We insert a complete set of coordinate bases into the one-body amplitude as 

(kle-^^-^lx) = Jdx a e~ lk - Xa (xale-^^-^lx) , (2-38) 

and adopt the semiclassical approximation in Eq. (I2-37|) on the right-hand side of 
the above relation. The classical action, S c (x a , t a ; x, to), is decomposed into two 
parts, 

Scip^ai tai "Kj ^o) — S c {x a) t a , x m , t m ) + iS c (x m , i m , a;, to) > (2 - 39) 

where x m (t m ) is a position on the classical trajectory. If we choose x m sufficiently far 
outside the range of the mean field potential, the first term is given by the classical 
action of the free motion, 

c / j. i \ fn{x a — x m ) 



and, substituting this form for Eq. (|2-38p and carrying out the integral with respect 
to x m , we obtain an expression for the amplitude, 

(k\e~ i ^( ta ~ t °^ \x) = j^ e ^ ik ' x ^ iE k{ t a-t m )+iS c {x m ,t m \x,t ) (2-40) 
If we assume that the amplitude, <Pk(to, x; t a ), is written in the form 

<Pk(to, t a ) = A k (t , x- t a ) e iS ^^) , (2-41) 
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a comparison to Eq. ()2-40p leads to a formula, 

S k (t , x; t a ) = -k ■ x - 5S c {x m , t m ; x, t ) - E k (t a - t ) (2-42) 

5S c (x m , t m ; x, t ) := S c {x m , t m ; x, t ) - {k ■ (x m - x) - E k (t m - t )} (2-43) 

and 

/ 2 

A k (t ,x;t a ) = y/2m(t - t m )/m A(x a ,t a ;x,t ) , 

where 5S c (x m , t m ; x, to) is the phase shift caused while a particle escapes from the 
range of the mean field potential. We note that the position x m (t m ) can be arbitrarily 
chosen on the classical trajectory, as long as t m is taken to be large so that the 
momentum becomes converged to the asymptotic momentum k. Our results do not 
depend on the dummy coordinate x m at t m . 

This procedure determines S k (to,x;t a ), if the set of boundary conditions, x and 
k, specifies a unique trajectory. In general, this uniqueness is, however, not guaran- 
teed mathematically, and there are some cases wherein more than two trajectories 
satisfy the same set of boundary conditions. In other words, the trajectories focus at 
an end of them, x. A similar difficulty is known in geometrical optics as the caustics. 
If we encounter this case, the standard semiclassical prescription breaks down, and a 
more elaborate treatment of the path integral is required to find a relevant solution 
of the amplitude. In our models examined in the following sections, the focusing of 
the trajectories does not arise in the source region, and the semiclassical approach 
works well. 

In the presence of the imaginary part of the potential, U, the prefactor may be 
given in the form 

lim Ak(to, x; t a ) = lim a(fc)exp 

t a — >CO t a — >0O 

This integral is calculated along the classical trajectory, x c (t), determined in the 
existence of the real part of the potential, V. A negative imaginary part provides 
the attenuation of the flux of particle while escaping the interacting region. 

Based on the amplitude given in Eq. (|2-4ip to Eq. (|2-44p . we first consider 
the modification of the one-body probability P\{k) in Eq. (|2T8p , Using the set of 
coordinates, x = {x\ + X2)/2 and r = x\ — X2, as in the preceding section, we obtain 

( -Pk(to,X 2 ]t a ) 

where the differential operator V acts on the coordinate x, and we have neglected the 
higher order derivatives of S k (to,x;t a ). This procedure is allowed for the potential 
that has a smooth spatial profile. This condition is satisfied automatically if the 
semiclassical approach itself is valid; the real part of the potential, V, should be a 
smooth function of the spatial coordinate over a wavelength of the particle. Owing 
to Eq. (|2-42|) . the first derivative of the classical action, X7S k (to,x;t a ), provides the 



ta 



dt U{x c (t)) 



/.i 



(2-44) 



= <p k (t , x + -r, t a ) * e^ r - VS ^ x ^ip k (t , x; t a ) (2-45) 
= cp k (t , x - X -v- t a ) ~ e^™* ' 3 ^ W*o, x; t a ) , (2-46) 
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initial momentum po at x (to) as 

VS k {t , x; t a ) = VS c (x m , t m ] x, t ) = p (*o, x, k) . (2-47) 

Inserting the above expressions, Eqs. (|2-45p and (|2-46p . into Eq. (|2T8p . we obtain 
the one-body probability given byEJ 

P 1 (k) = JdxJdp5(p- Po )f(to,x,p)\A\ 2 (2-48) 

= \a{k)\ 2 jd A x S{x, Po ) e - 2 ^' x 'V , (2-49) 

where S(x,p) is the source function defined in Eq. ()2-23j) . and "f(to, x,k) = — J^dtU{x c {t)) 
provides the amount of the absorption. For the sake of simplicity, we have suppressed 
the arguments of the initial momentum po(to, x, k) given in Eq. (|2-47p . 

This result agrees with our expectation: the real part of the potential causes a 
shift of the momentum, and the imaginary part provides an attenuation factor of 
a particle traversing the cloud. Note that the shift of the momentum distribution 
is a classical effect that is also obtained from the classical treatment GEJ'GU j n the 
quantum description, the phase shift of the amplitude caused in the interacting region 
contains the information on the momentum shift. The absorption originates in the 
deviation from the classical trajectory due to the collision in the cloud. This effect 
is mainly attributed to the formation of p meson resonance in the pion cloud, and 
may be partly due to the formation of A and K* with much rarer nucleon and kaon, 
respectively. In HBT interferometry, the modification of the one-body probabilities 
reflects in the normalization of the correlation function, C. 

We proceed to the joint probability P2(k\, k 2 ), and examine how the interference 
pattern is distorted owing to the modification of the amplitude. We begin with 
the expression given in Eq. (|2-19p . which is obtained by the factorization of the 
amplitude and density matrix of the particle pair preserving the boson symmetry. 
In this relation, we have the interference term, 

F(fci,fc 2 ;<a) = y *ci /" iaB2 <*fi I iy"(t a , to) |aci> <asi|p(*o) |aca> <asa 1 1?" 1 " *o) |As 2 )(2.50) 

in which the effects of the interaction are incorporated in the amplitudes examined 
above. Inserting the amplitudes given in Eqs. (|2-45|) and (|2-46|) into Eq. (|2-50|) . 
we find that the interference term is expressed with the difference in the classical 
actions 

F(k 1 ,k 2 ;t a ) = a(k 1 )a(k 2 ) JdxJdpS (p - p ) f(t ,x,p) 

x e -7(*o>a5j*l)-7(*o,a!,fe2) gif^fej (to,x;t a )-S k2 (t ,x;t a )} (2-51) 

where the initial momentum po(to,x,k) is given by the derivative of the classical 
action as in Eq. (|2-47p . Here, we expand the actions with respect to the relative 
momentum q, assuming that the interference arises in the small q regime compared 
with the magnitude of the averaged momentum k. This approximation is valid for 
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the interferometry in the ultrarelativistic heavy ion collision, in which the averaged 
momentum is in hundreds of MeV, while the relative momentum is on the order of 
the inverse of the source size ~ 10 fm, which is in tens of MeV. Up to the leading 
order of the relative momentum, we have 

<Sfei(*o,sc;*a) - Sk 2 (t ,x;t a ) ~ q ■ V k S k (t , x; t a ) , 

where V& operates on the averaged momentum k explicitly, and the energy E k 
implicitly. The attenuation factor 7 is also expanded as j(to, x , &i) + 7(io> x , k 2 ) — 
2'y(to,x,k) . Using these expressions, we obtain 

F{kx,k 2 ;t a ) = a{kx)a(k 2 ) Jdx f(x,p ) e - 2 ^ to ' x ' k) e ^^ k {t ^t a ) (2 . 52) 

Note that the emission occurs at the time to(x) depending on the position x. 

To see how the source images are distorted by the mean field interaction, we 
write 

S k (t , x; t a ) = S ok (t , x; t a ) + 5S k (t , x) , 

where So k (to, x; t a ) is the classical action of the free motion given by Eq. (|2 42|) as 
S'ofe(to) x; t a ) = —k ■ x — E k {t a — to), and 5S k (to, x) is the phase shift caused by the 
mean field interaction, 

5S k (t ,x) = S k (t , x; t a ) - S ok (t , x; t a ) . 

The difference in the actions, or the phase shifts, does not depend on the time t a , 
since the difference is generated only in the interaction region localized in the vicinity 
of the source. Taking the derivative with respect to k, we find a relation 

VfeSfc^o, x; t a ) = -x - v k (t a - t ) + V k 5S k (t x) , (2-53) 

where, as in the previous section, the velocity is defined by the off-shell momentum, 
v k = V k E k . We, however, usually assume the on-shell condition, E\ k 2 + m 2 , 
for the small relative momentum. By inserting Eq. (|2-53p into Eq. (|2-52p . the 
interference term can be expressed as 

F(k 1 ,k 2 ;t a ) = a(k 1 )a(k 2 ) 

x jdx f(x, Po ) e - 2 7(*o,=,fe) e -iq-(x-v k t -V k 5S h (to,x)) _ (2-54) 

Note that the derivative of the action, V k 5S k (to, x), is not the variation of the action 
with respect to the momentum shift on a classical trajectory, but originates in the 
difference in the two actions of the adjacent trajectories, labeled by the asymptotic 
momenta, k\ and k 2 . 

Together with the one-body probability in Eq. (|2-49p . the correlation function 
defined in Eq. (|2-1|) is then given by 



C=l + 



Jdx f(t ,x,p ) e - 2 T(*°,*,fc) e -fc>-«kto-Vk*5*(to,«>)) 2 



Jdx f(t ,x,p ) e- 2 T(*°- a! ' fe ) 



(2-55) 
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or, using the relation (|2-27p . we obtain 

Jd 4 x S(x,po(x,k)) e- 2 ^' fc ) e iq^+d»ss k {x)) 



C = 1 



Jd 4 x S(x,p {x,k)) e - 2 7(^) 



(2-56) 



where x is the space-time coordinate, x M = (to,x). The source function S(x,p) is 
defined in Eq. (|2-23p . In these expressions, we have approximated the momenta 
appearing in the denominator as k\ ~ k 2 — k for the small relative momentum. The 
temporal component of the derivative d£ operates on the energy E k , and the spatial 
components operate only on the explicit momentum dependence k, but not on the 

energy E k . 

In Eq. (|2-55p . or Eq. ()2-56f) . we find the two modifications by the real part 
of the mean field potential: the momentum shift in the source function due to the 
acceleration, and the coordinate shift in the kernel of the transformation given by 
the derivative of the classical action. As in the one-body probability, the imaginary 
part provides the attenuation factor. Incorporating these effects and assuming the 
Gaussian shape of the source profile as in Eq. (|2-32|) . we still find the set of formulae 
in the form given by Eq. (|2-33|) . The modifications merely reflect in the definitions 
of the deviation x M and the expectation value (&), 

x» = ( X M + d£6S k (x)) - (x» + d%6S k (x)) (2-57) 
_ f0S(x,p o ) e'^^d'x 

{ } ~ JS(X, P0 ) e -27(*,fc) d 4 x ' ^° 6 > 

in which we find the three effects mentioned above: the momentum shift, the coordi- 
nate shift and the attenuation. The arguments of the initial momentum Po(x, k) are 
suppressed in Eq. (|2-58|) . Replacing Eqs. (|2-34p and (|2-35|) by the above definitions, 
we obtain the Gaussian parameters in the presence of the mean field interaction. 

To examine the distortion of the HBT image more intuitively, we define the 
apparent emission point x' as 

x' = x - v k t (x) - X7 k 5S k (t (x),x) , (2-59) 

where the right-hand side of the above relation appears in the kernel of the integral 
transform in Eq. (|2-55p . Transforming the integral variable x to x', we obtain the 
integral for the interference term given by 

F(k q ,k 2 ;t a ) = a(fc 1 )a(fc 2 )e-^ ta fdx 1 f eS (x',k) e^ x ' (2-60) 

with the definition, 

f eS (x',k) = J(x,x';k) f(t ,x,p ) e - 2 ^'.^) , (2-61) 

where J(x,x';k) = d(x, k)/d(x', k) = [d(x',k)/d(x,k)]~ 1 is the Jacobian of the 
transformation in Eq. (|2-59j) . Note that the original coordinate x appearing in the 
expression of f e g(x',k) is understood to be a function of x' and k via Eq. (|2-59p . 
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We have obtained the static case of the above relations in the previous workp^ in 
which the temporal term Vkto(x) in Eq. (|2-59p is absent. The modification of the 
propagation after the emission, arising as the phase shift of the one-body amplitude, 
is transferred to a distortion of the spatial distribution by the apparent shift of the 
emission point defined in Eq. (|2-59p . We find that the correlation function is written 
in the same form as in the free case using the Fourier transform, 



C(k,q) = l + 



Jdx'p eS (x',k)e i ' 1 - 



(2-62) 



with the effective distribution function defined by 



p«&> k ) = r . f( / cG{x, : k) 2y(t0Xk) ■ ( 2 -63) 

jdx j(to,x,p ) e -Wo,*!*) 

Owing to the nonlinear nature of the transformation in Eq. (|2-59p , the original coor- 
dinate x is mapped onto the apparent coordinate x' with a deformation of the grids, 
and it induces the nonlinear mapping of the source distribution function onto the 
effective one defined on the apparent coordinate system. This nonlinear nature acts 
to distort the image viewed on the coordinate x'. The attenuation factor provides 
the opacity of the cloud. Integrating both sides of Eq. (|2-63[) with respect to x/, we 
observe that the effective distribution function p e g(x' , k) is normalized to be one at 
each k even in the presence of the absorption, whereas its effect reduces f e g(x', k) in 
Eq. (|2-6ip . We find that the normalization of the effective distribution function is 
compensated by the single-particle spectrum in the denominator of Eq. (|2-63p . which 
contains the same attenuation factor. The effective distribution function represents 
the distribution of the emission points of the surviving particles. The absorptive ef- 
fect does not change the magnitude of the correlation term at the vanishing relative 
momentum, C(k,0) = 2. In we illustrate the distorted profiles of the effective 
distribution function using a schematic model of the potential. 

We briefly comment on the computation procedure concerning the following two 
points. In Eq. (|2-59|) . we take the derivative of the classical action with respect to 
k using the cylindrical coordinate in the momentum space, since the conventionally 
used coordinate refers to the direction of k. The apparent shift of the emission point 
is represented by its components as 

x > = x - v±to (x)-™%&, 

y' = y-\ 9 -W^ (2-64) 

, _ _ dss k (£) 

- Z ~ Z dz ' 

where k, 0^ and v± are the magnitude and angle of momentum k, and the magnitude 
of the velocity Vk, respectively. Owing to the definition of the coordinate system in 
which 

Vy,Vz) = (v±, 0, 0), the temporal term proportional to the velocity 
is preserved only in the outward coordinate x. Even in the free case without the 
derivative terms, we find the apparent extension in the outward direction due to the 
presence of the temporal term. This effect is absent for the static source and the 
instantaneous emission, which do not have the temporal structure. 
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To calculate the classical action of the trajectory specified by the emission point 
x and the asymptotic momentum k, we have to find the angle of the initial mo- 
mentum p(t = 0) such that, it satisfies the matching condition at the large time 

as 

p(oo) = k . 

Owing to the energy conservation law, the magnitude of the initial momentum p(0) is 
simply related to that of the asymptotic momentum k. On the other hand, the angu- 
lar momentum conservation is not adequate to fix the angle of the initial momentum, 
and we obtain the angle by explicitly solving the equation of motion. Regarding the 
left-hand side of the above condition as a function of x and p(0), we solve it with 
respect to the angle of p(0). A good number of the trials of computing the initial 
value problem, and their interpolation, enable us to find the solution. Carrying out 
such calculations for the adjacent values of the momentum k, we obtain the deriva- 
tive of the phase shift, bx = — V^dS^x, k), and the apparent shift of the emission 
point defined by Eq. (|2-59p . 

§3. Profiles of the distorted images 1 ' 1 

In this section, we examine the distortion of the HBT images based on the 
effective source distribution function defined in Eq. (|2-63|) for the static source. Here, 
the temporal term in Eq. (|2-64p is absent. We show how the difference in the effects 
of the mean field interactions, i.e., attraction and repulsion, appears in the distorted 
images. Nonrelativistic analysis is performed using a schematic model of the mean 
field potential given below. Recall that the coordinate system on the transverse 
plane is spanned by the direction of the averaged momentum of particle pair k and 
the direction perpendicular to it. They are called outward (x) and sideward (y), 
respectively. Roughly, the particle pair is detected in the positive direction of the 
outward. 

For simplicity, we assume the factorization of the bare source distribution func- 
tion as 

f(x,p)= p(x)f th (p;T,p) , (3-1) 

where p(x) is the spatial profile of the source and f t h(p',T, p) is the thermal spec- 
trum specified by temperature T and pion chemical potential p. We input a two- 
dimensional Gaussian profile on the transverse plane, p(x) = po e~ x / 2A , with A = 5 
fm and the Bose-Einstein distribution function with T = 140 MeV and p = MeV. 

We have chosen a schematic pion optical potential assuming a two-range Gaus- 
sian shape for the real part: 

V(x) = Vie + V 2 e 2A 2 , (3-2) 
and another Gaussian shape for the imaginary part: 

U(x) = - U e 2A o . (3-3) 



18 



K. Hattori 
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Fig. 1. Original grids mapped onto the apparent coordinate system: Curves on the left (right) 
panel indicate the deformation of the grids by repulsive (attractive) mean field interaction. We 
show the deformation at the asymptotic momentum k = 100 MeV in the range \x\ < 10 fm and 
\y\ < 14 fm. The deformation induces the apparent distortions of the source images shown in 
Figs. Hand [3] 



The shorter ranges, Ai and Ao, are taken to be consistent with the source profile A 
given above. The other range A2 gives the thickness of the meson cloud. We have 
chosen A x = 5 fm, A 2 = 5 fm and a = 10 fm; V% = ±10 MeV, V 2 = ±2 MeV and 
Uq = 0.1 MeV. We compare the effects of the repulsive interaction and the attractive 
one. 

In £ )2,2l we have shown that the phase shift caused by the real part of the mean 
field potential results in the apparent shift of the space-time coordinate. We have 
defined a nonlinear transformation in Eq. (|2-64|) , which maps the original coordinate 
system onto the apparent one. In Fig. [H we show the location of the original grids 
mapped on the apparent coordinate system at the asymptotic momentum k = 100 
MeV, and in the range \x\ < 10 fm and \y\ < 14 fm. Curves on the left (right) 
panel represent the original Cartesian coordinate under the effect of the repulsive 
(attractive) mean field interaction. Induced by the deformation of the grids, the 
profiles of the image are distorted on the apparent coordinate system, as shown in 
Figs. [2] and [3j The deformation of the grids in Fig. [1] represents the geometrical 
interpretation of the effects by the real part of the potential. However, note that 
each curve does not indicate a refracted trajectory of the emitted particles, which is 
often described to express the lensing effect in geometrical optics. 

Using the contour plot, we show the profiles of the effective distribution function 
defined in Eq. (l2 T 63l) at k = 100 MeV in Fig. El and k = 150 MeV in Fig. GJ The 
upper panels in each figure show the distortion of images without the absorption, 
whereas the lower panels contain the effect of the absorption with IAq = 0.1 MeV as 
well as the effect of the real part of the potential. The two panels in the left column 
show the cases with the repulsive interaction, and those in the right column the 
attractive one. The vertical and horizontal axes extend into the apparent outward 
and sideward directions, respectively: thus, the location on each panel is the apparent 
emission point. In Figs. [5] and [31 the distorted perimeters of the plot range indicate 
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Fig. 2. Contour plots of the effective source distribution function, p c g(x' , k), defined in Eq. (|2-63p . 
The averaged momentum of the pair is chosen to be k = 100 MeV. The panels in the left (right) 
column show the distorted profiles due to the repulsive (attractive) interaction. We incorporate 
the effect of the absorption in the lower two panels. The vertical and horizontal axes are the 
transformed outward and sideward coordinates, respectively. We find the deformation of the 
grids by the nonlinear coordinate transform in the distorted perimeters of the plot range. 



the deformation of the original grids. 

We first compare the cases without the absorption shown in the upper row. At 
k = 100 MeV, we find that the apparent source image is elongated in the outward 
direction by the repulsive interaction, while it is stretched in the sideward direction 
by the attractive one. These effects are pronounced in the backside region of the 
source, since the amplitudes of the particles, penetrating into the interior of the 
source, are strongly distorted. However, the absorption damps the amplitudes of 
such particles emitted in the deep interior and backside region of the source, and 
diminishes the remarkable effects caused by the real part of the mean field potential. 
On the lower two panels, we find the images stretched in the sideward direction 
with the dominant distributions on the detector side of the source. The cylindrical 
symmetry assumed for the original source profile is broken by the interaction, since 
the strengths of these effects depend on the emission point. As we found in Figs. 
and [3l only the reflection symmetry with respect to the outward axis survives in the 
presence of the mean field interaction. At k = 150 MeV, these effects are found in the 
same trend, but are less effective: the distortion of the amplitude is relatively weak 
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Fig. 3. Contour plots of the effective source distribution function, p e g(x',k), at k = 150 MeV. In 
each panel, we incorporate the same interaction as in the panel at the same location in Fig. [5] 



at the high momentum, if the mean field potential does not depend on momentum. 

In our analysis, the origin of the distortion, by the real part of the mean field 
potential, is the nonlinear mapping of the coordinate system. The transformation is 
given by the derivative of the phase shift with respect to the averaged asymptotic 
momentum, k = (k\ + k2) /2. Recall that the derivative indicates not the variation of 
the action with respect to the momentum shift on a trajectory, but the difference in 
the actions of two adjacent trajectories. In the presence of the one-body mean field 
potential, we have different phase shifts on each trajectory labeled by fci and k2, 
and the common emission point x. Accordingly, their relative momentum is shifted 
by the one-body potential due to the coordinate dependence of the force to which 
each particle is subjected. The shift of the relative momentum is the origin of the 
distortion of the images observed utilizing a quantum interference. Thus, we remark 
that the distortion of the images in HBT interferometry is not a direct analogue of 
that in geometrical optics described on the basis of refracted trajectories. 

To emphasize an aspect of the quantum interference, we examine the phase 
shift using a Glauber-type approximation without the refraction of trajectories. It 
assumes the straight line trajectory in the interaction region, so that the phase shift 



Distortion of HBT Images 



21 





-20 -rs -ro 



to rs 20 




-20 -rs -ro 



ro rs 20 



Fig. 4. Effective source distribution function, p e g (as', k), at k = 100 MeV obtained with the 
Glauber-type approximation: We have chosen the same potentials as those in Fig. These 
results qualitatively reproduce the features found in Fig. [2] 



is given by a simple formula: 

m MT) 

5S Glanher (x,k) ~ - — / V(Vu 2 + 6 2 ) , (3-4) 

i«(o) 

where 6 = \Jx 2 — (x • k) 2 /k 2 is the "impact parameter" of the trajectory defined 
by the distance between the straight line and the center of the potential. Position, 
u(t) = x(t)-k/\k\, is the coordinate along the straight line at time, t. We then obtain 
the coordinate transformation (|2-64p in the outward direction from the dependence 
on the magnitude of momentum, and in the sideward direction from the dependence 
on argument 6k through b. The distortion of images at k = 100 MeV is shown in 
Fig. HI using the same potentials as in Fig. [2j We find the qualitatively same results 
assuming the straight line trajectory. 

§4. Mean field interaction in the freeze-out stage 

4.1. Phenomenological mean field interaction 

In this section, we examine the dispersion relation of the pion in the cloud formed 
by the other evaporating particles. In analogy to the optical model employed in the 
analysis of nuclear reactions, we construct a phenomenological self-energy based on 
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the two-body forward scattering amplitude. Although an escaping pion possibly 
interacts with other pions, nucleons, kaons and resonances directly emitted on the 
freeze-out hypersurface, we focus on the interaction with the most abundant pion in 
this work. The magnitude of the mean field interaction and the sign of its real part 
are determined using the results of a phenomenological analysis on the two-body 
7T7T scattering amplitude,®'® which fits the measured amplitude incorporating the 
constraint of the underlying chiral symmetry. The presence of the p meson resonance 
peak plays important roles in both the real and imaginary parts of the mean field 
interaction. 

The modification of the dispersion relation is expressed with the self-energy 77: 

E 2 p = p 2 + m 2 + 77(p) , (4-1) 

where p and m are the momentum and mass of the pion escaping from the cloud, 
respectively. Since the system becomes diluted in the freeze-out stage owing to the 
evaporation, two-body scattering of thermal pions dominates the contribution to 
the pion self-energy. In Fig. [SJ we sketch the scattering with a thermal excitation, 
where the closed line indicates an excited medium pion carrying momentum p'. Up 
to the two-body scattering of on-shell excitations, the diagram contains only a single 
thermal loop, and the gray blob represents the forward scattering amplitude in the 
vacuum. The pion self-energy 77 described in Fig. [5] is written as 

n h (p) = - £ [-^2n5(p' 2 -m 2 )T^ 2 (s) f l2 (p') , (4-2) 
r a =o,±i ^ ' 

where the isospin component of the medium pion, 72, is summed over -1,0 and +1, 
and s is the Mandelstam variable, that is, the squared center-of-mass energy of the 
two-body scattering defined by s = (p + p') 2 . T Ill2 (s) and f l2 (p') are the two-body 
forward scattering amplitude in the vacuum and the distribution function of the 
medium particle carrying the momentum p' and the isospin 7 2 , respectively. The 
overall sign is up to the definition of the phase shift in T Ill2 (s); we take the minus 
sign so that we have a negative phase shift for repulsive scattering. The distribution 
function and the on-shell condition come from the thermal part of the pion propa- 
gator. Here, we suppose that the renormalization has already been achieved in Eq. 
(|4-2|) . and input the physical pion mass, m = 139.5 MeV, and the phenomenological 
amplitude into Eq. ()4Tp . The integral with respect to the momentum p' does not 
pose further ultraviolet divergence owing to the distribution function. 

Neglecting the small difference in the light quark masses, the scattering ampli- 
tude by the strong interaction is blind to the third component of the total isospin. 
Here, we assume that each species of pion is produced in the same amount in the 
RHIC experiment GSIES which we write as f(p') = f l2 {p') for 7 2 = 0,±1. In the 
isospin symmetric limit, the right-hand side of Eq. (|4-2p is then simplified to 

nip) = -J-0f 2 ^p' 2 ~ ™ 2 ) T (s) M , (4-3) 
too = 3 Qt ( s ) + |z\0O + jjr 2 ( fl )) , (4-4) 
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Fig. 5. Mean field interaction with thermal exci- Fig. 6. A model profile of the freeze-out proper 
tations: In the pion cloud, an escaping pion time: On the transverse plane, the "freeze- 
picks up a thermal pion carrying the momen- out" in the surface region precedes that in the 
turn p' and the isospin 12- The gray blob interior. Thus, the freeze-out surface shrinks 
represents the forward scattering amplitude in as time advances, 
the vacuum. 



where, owing to the isospin symmetry, the self-energy is independent of the isospin 
carried by an escaping pion. The sum over the isospin component of the medium 
pion in Eq. (|4-2p results in the overall multiplication factor attached to the isospin- 
averaged scattering amplitude T{s) in Eq. (|4-4p . The subscripts on the right-hand 
side of Eq. (|4-4p denote the total isospin channels. To evaluate the self-energy given 
by Eq. (|4-3p . we need two distinct physical quantities, T(s) and f(p'), which provide 
a microscopic foundation to the mean field interaction. 

We first examine a simple time evolution of the distribution function f(p'). The 
dynamics of the expanding source is important in two respects. First, the density of 
pions, strolling in the vicinity of the source, becomes diluted owing to the evapora- 
tion, which diminishes the modification of the dispersion relation as time advances. 
In the time evolution of the matter created in the ultrarelativistic heavy ion colli- 
sion, the evaporation is promoted by the transverse expansion and the approximately 
boost invariant strong longitudinal expansion. Second, the collectivity changes the 
typical scattering energy of pions compared with the case without it. This effect re- 
flects in the self-energy through the Mandelstam variable in the scattering amplitude, 
T(s). 

To reflect the space-time geometry of the source, we adopt a phenomenological 
extension to Eq. (|4-3p . The distribution function is here assumed to depend on 
the space-time coordinate as well as the momentum. Assuming the free streaming of 
medium pions for the first approximation, the self-energy is then obtained by solving 
the free classical transport equation: 

n(t,x,p) = -J0z 2v5{p' 2 -m 2 ) T(s) f(t,x,p') (4-5) 

p'^f(t,x,p') = S(t,x,p') (4-6) 

where the self-energy II(t, x,p) depends on the momentum p through the center-of- 
mass energy, y/s. A source term S(t,x,p') describes the emission at the freeze-out 
hypersurface. Solving Eq. (|4-6p under the initial condition, f(t = 0,x,p') = 0, we 
obtain the distribution of the pion after the emission. 
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To discuss the longitudinal boost invariance, it is convenient to introduce the 
kinematical variables defined by 



m ± = p 2 ± + m 2 , \ t = Vt 2 - z 2 

where pj_ is the magnitude of the transverse momentum, p± = (p x ,Py)- Using these 
variables, we rewrite Eq. (|4-6|) to exhibit the manifest symmetries as 

9 1 , d 1 _ \ „. S(t, x,p) /A „, 

tanh£— + — p± -V± ) f(T,r,p ± ,u,t )= V , 4-7 

or t c*4 mj_ cosh £ / mj_ cosh £ 

where we dropped the primes on the momenta. We have the two differences £ = Y— r/ 
and uj = (f) — 6, where the angles, and 0, are the arguments of the transverse 
coordinate x± and the transverse momentum respectively. These differences are 
invariant under the Lorentz boost in the longitudinal direction and the rotation on 
the transverse plane, respectively. We find the manifest cylindrical symmetry on 
the left-hand side of Eq. (|4-7p . noting that the inner product of the vectors on the 
transverse plane is given by the difference in the angles, p± • Vj_/ = |p_i_||Vj_/| cos a;. 
Thus, the cylindrical symmetry of the distribution function is preserved under the 
time evolution, as long as the source function has the same symmetry in the form, 
S(t, r,p±,u,g). 

Provided that the evaporating pions have the thermal spectrum on the freeze-out 
hypersurface, the source functions S(t, x,p) in Eqs. (|4-7p and (|2-23[) are consistently 
given by 

S{t,x,p) = — Lr ! fe q (x',p) 5\x-x') p v da^x>) , (4-8) 



(2tt) 3 J JeqK 
1 

exp{{p v u 1/ - n)/T} - 1 



where x' is the emission point on the hypersurface, and da^x') is the normal vector 
of it. The integral is carried out over the hypersurface owing to the delta function. 
The thermal distribution function f eq (x,p) is specified by the macro variables: tem- 
perature T, pion chemical potential \x and the flow vector it^. The longitudinally 
boost invariant expansion is described using the Bjorken flowP^ 

u ^ x ) = —^=( t -, v ^ Z -\ , (4-10) 
'l-vl \ T T ) 

where v± is the transverse flow vector. Its magnitude depends on the transverse 
radial coordinate, r. The energy on the rest frame of the fluid cell is written with 
the flow vector as 

p»u v = m±COShC -^ COSa; , (4-11) 
-v 2 
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where p±_ and v± are the magnitudes of the transverse momentum and the transverse 
flow velocity. Owing to the above expression, we find the boost invariance and 
rotational invariance of the source function in Eqs. (|4-8|) and (|4-9|) . The normal 
vector of the hyper surf ace is given by 

/ dtf(r,z) \ 

da^{x) = Tf I cosh rj, cosh 77 — ^— - — e_|_,sinhr/ I rdrdOdi] , (4-12) 



where the freeze-out time, tf(r,z) = y T f( r ) +z 2 , is specified by the r-dependent 

freeze-out proper time Tf(r) given below. On the transverse plane, the "freeze-out" 
may proceed from the surface region to the deep interior due to the finite source 
volume and the finite sound velocity. 

Motivated by the hydrodynamical simulations DSiED we have used the following 
parameters and analytic functions in the numerical analysis. We take T = 130 MeV 
and \i = 30 MeV for the temperature and the pion chemical potential, respectively. 
The profile of the radial flow is approximated using a linear function of the radial 
coordinate r, \v±\ = 0.06 r c , and the freeze-out proper time using an ellipsoidal 
curve, Tj(r) = Tiif e yl — r 2 /rg. The two radii, riife and ro, are roughly the life- 
time of the source and the radius of the cylinder soon after the nucleus collision, 
respectively(see Fig. [6]). We have chosen riif e = 15 fm/c and ro = 8 fm. 

Concerning the classical transport equation in Eq. (|4-6p , we should comment on 
the following two points. First, the source term introduced to describe the emission 
on the freeze-out hypersurface can be negative at the timelike part of the hypersur- 
face, owing to the longstanding issue of the Cooper-Frye formula. Its contribution 
induces a negative value of the distribution function f(t,x,p). To cope with this 
issue, we have confirmed to what extent it could reflect in our results. If we find the 
negative distribution, we cut off its contribution to the integral in Eq. (|4-5p . Com- 
paring the results obtained with and without the negative distribution, we observed 
that the self-energy changes within 2%, and the final results, the Gaussian radius 
parameters in Fig. [13] and the transverse spectrum in Fig. [T4"l change by not more 
than 1%. We consider that the integral in Eq. f|4-5j) with respect to the configura- 
tion of medium pion smears the negative contribution. The negative contribution 
is negligibly small, at least within our choice of the parameters and the profile of 
the freeze-out hypersurface specified above. Although the longstanding issue of the 
Cooper-Frye formula still remains, our results are barely affected. 

Second, to treat the motion of medium pions self-consistent to the effects of the 
mean field interaction, it will be necessary to incorporate the Vlasov term and the 
collision term into the transport equation in Eq. (|4-6p . They are described using the 
self-energy II(t,x,p), which couples the transport equation (|4-6p to the expression 
of the self-energy in Eq. (|4-5p . However, owing to the following discussions, we again 
deduce that the small deviations in the configuration of medium pions are smeared in 
the integral in Eq. (|4-5|) . as is observed for the contribution of negative distribution. 

The Vlasov term describes the variation of the distribution due to the momentum 
change, or the deviation of the classical trajectory from the free motion, by the 
external force. We find in the equation of motion (|5-7|) that the force is proportional 
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to the gradient of the real part of the self-energy. The validity to neglecting the 
Vlasov term is not obvious a priori. Nevertheless, in the computation of the single- 
pion spectrum shown in Fig. [T~H we find that the shift of the transverse momentum 
by the mean field interaction is less than 20 MeV. Provided that the momentum shift 
is caused in the time scale on the order of the duration, rnf e , the deviation of the 
classical trajectory from the free motion is estimated to be \Ax\ ~ Tife ^ 1 fm. 
The momentum shift and the deviation of the trajectory are almost an order smaller 
than the magnitude of the momentum and the extension of the source region, tq, 
respectively. 

The collision term describes the loss and gain of the pions in the phase space 
volume element. The optical theorem relates the imaginary part of the forward 
scattering amplitude T(s) to the total cross section of the binary collision, tttt — > 
{anything}. Thus, the loss of the pions distributed in the interval, x ~ x + Ax 
and p ~ p + Ap, at t is described by the imaginary part of the self-energy as 
Im [n(t,x,p)] f(t,x,p) up to the normalization. In Fig. [7J we find that the imag- 
inary part of T(s) is dominated by p meson resonance peak. Therefore, the loss is 
mainly attributed to the formation of p meson resonance, 7r(p)7r(j/) — > p(p + pi) F^l 
The gain is dominated by the inverse reaction, that is, the decay of p meson reso- 
nance p{p + pi) — > tt(pi)tt(p2), which has the almost 100% branching ratio. Because 
this reaction, n(p)ir(p') — > p(p + p/) — > vr(pi)7r(p2)> appears in the p-wave scatter- 
ing, the differential cross section of this process is proportional to the square of the 
Legendre function, |Pf = i(cos 6)\ 2 = cos 2 9, where 9 is the scattering angle in the 
center-of-mass frame. The momenta in the final state are localized in the forward 
direction: therefore, we have p ~ p\ and p' ~ P2, or p ~ p2 and p' ~ p±, owing to 
the momentum conservation. Although the finite width of the angular distribution 
in the final state causes a small variation of the distribution function, this effect is 
smeared in the integral in Eq. (|4-5|) . As far as we use the distribution of the free 
streaming pions to obtain the self-energy, we expect that the effects of the Vlasov 
term and collision term are, in practice, not relevant to the purpose of this work. 
However, a self-consistent framework with these terms will be necessary to describe 
possibly more efficient effects of the mean field interaction in the denser system. 

Next, we show the momentum dependence of the forward tttt scattering am- 
plitude in the vacuum. Pion is known as an approximate Nambu-Goldstone boson, 
which has been studied on the basis of the underlying chiral symmetry, and its break- 
ing. Here, we use a phenomenological analysis of the tttt scattering amplitude in the 
s-wave and p-wave. It was achieved on the basis of the Roy equation™ which is the 
manifestation of the fundamental properties of the scattering amplitude: analyticity, 
unitarity and crossing symmetry. This analysis incorporates the constraint by the 
chiral symmetry, and precisely reproduces the behavior of the p meson resonance 
peak. 

In each total isospin channel, the scattering amplitude T 7 (s,i) is decomposed 



*' The absorption by the mean field interaction is thus microscopically originated in this reac- 
tion. 
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Fig. 7. Elastic forward irir scattering amplitude T I (s,0) in the s-wave and p-wave: The panels 
on the left and right show the real and imaginary parts of the amplitude, respectively. The 
horizontal axes show the center-of-mass energy, yfs, of the two-body scattering. Curves show 
the amplitudes in each isospin channel, T°, T 1 and T 2 , and their taken with the degeneracy 
factor for the isospin multiplet, as shown in Eq. (|4-4[) . We find a narrow peak in the imaginary 
part of the 7=1 channel around y/s — 770 MeV, located on the mass of p meson resonance. 
The experimental data^ are indicated by markers for the isospin channels. 

into the partial waves as 



where the Mandelstam variables, s and t, are related to the energy and the scattering 
angle in the center-of-mass frame. Below the threshold of two kaons, ^/s nn ~ 1 GeV, 
we focus on the elastic scatterings without the inelasticity. The analysis provides 
the phase shift SUs) in the s-wave and p-wave. 

We show the forward scattering amplitude T (s, 0) in Fig. [JJ The panels on the 
left and right show the real and imaginary parts of the amplitude, respectively. The 
horizontal axes show the center-of-mass energy of the two-body scattering. Curves 
display the amplitudes in each isospin channel, T°, T 1 and T 2 , and their average 
taken with the degeneracy factor for the isospin multiplet as shown in Eq. (|4-4j) . 
The maximum value in the 1 = 1 channel is three times larger than those in the 
other two channels owing to the degeneracy factor for the angular momentum, t. 

In the s-wave, the real part of the amplitude in the 1 = channel has an 
attractive regime at low energy, while that in the 1 = 2 channel exhibits weak 
repulsion. The phase shift passes ir/2 at the intercept of the horizontal axis on the 
left panel, or the center of the peak on the right panel. The amplitude in the 1 = 
channel shows a slow variation against the increase in y/s, and the imaginary part of 
it has a broaden peak. The energy dependence in the 1 = 2 channel is much weaker, 
and the phase shift even does not pass ir/2 in the range shown in these plots. On 
the other hand, the amplitude in the 1 = 1 channel has a strong energy dependence. 




(4-13) 



(4-14) 



(4-15) 
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t=5(fm/c) t=12 (fm/c) t= 19 (f m/c) 



Fig. 8. Time evolution of the effective pion mass defined by the real part of the self-energy in Eq. 
(|4T6[I : We show the two-dimensional spatial profiles of the effective pion mass at the transverse 
momentum, p± = 150 MeV. In this regime, the mean field interaction acts on the HBT images 
most effectively. We find the effects of the attractive mean field interaction in the vicinity of the 
freeze-out hypersurface, in the exterior of which any interactions are assumed to be absent in 
the hydrodynamical modeling. The sharp boundary is smeared by the mean field interaction, 
and the smeared boundary shrinks as time advances. We do not find the effect of the repulsive 
mean field interaction in these plots. 




p =100 (MeV) p =300 (MeV) p. =500 (MeV) 



Fig. 9. Momentum dependence of the effective pion mass: These panels exhibit the transverse 
momentum dependence of the effective pion mass at the moment, t = 5 fm/c. At p± = 500 MeV, 
we find effect of the repulsive mean field interaction in the backside region of the source, where 
the center-of-mass energy of irn scattering is typically large owing to the collective expansion of 
the medium. 

Owing to the optical theorem, the narrow peak in the imaginary part is related to 
the peak in the cross section, and its location at y/s ~ 770 MeV is identified to the 
mass of p meson resonance. While the average in the s-wave scattering displayed 
with a thin solid line exhibits a weak repulsion, the contribution of p-wave scattering 
overwhelms it to provide an attraction below the mass of p meson, as indicated by 
the thick solid line. On the right panel, the thick solid line for the average is close 
to the dash-dotted line for the 1 = 1 channel. The absorptive effect described by 
the imaginary part is dominated by the formation of the p meson resonance in the 
1 = 1 channel. 

Carrying out the convolutional integral in Eq. (|4-5p . we obtain the phenomeno- 
logical self-energy. To illustrate the effect of the real part of the self-energy, we show 
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Fig. 10. Contour plot ol the attenuation factor for the pion emitted at the position (x,y) and 
detected with the asymptotic momentum The attenuation becomes stronger as the color 
becomes darker. As given by Eq. (|5T0f) . the integral of the imaginary part of the self-energy 
is carried out along the world line, which is computed under the effect of the real part. Each 
position on the panel is the emission point located on the freeze-out hypersurface viewed on the 
transverse plane. Note that the emission at each time is plotted on the same panel. 



the effective pion mass, which is defined by 



ill. 



■es(t, x, p) = yjm 2 + Re77(t, x,p) . (4-16) 



In Figs. [5] and O we show the spatial profiles of the effective pion mass with its 
time dependence and momentum dependence, respectively. Note that we focus on 
the central rapidity, p z = 0, and that the argument of the transverse momentum p± 
carried by the escaping pion is fixed at <j> = without loss of generality 

In Fig. we show the time dependence of the effects of the mean field interac- 
tion, where time advances from left to right, and the transverse momentum is fixed 
at p±_ = 150 MeV. We find the reduction of the effective pion mass in the vicinity of 
the freeze-out hypersurface at each moment. Note that, since we do not consider the 
interactions in the interior of the hypersurface, we have a flat region around the ori- 
gin with m c fj ~ 139 MeV, which shrinks as the "freeze-out" proceeds. The reduction 
of mass is the effect of the attractive mean field interaction obtained on the basis 
of the microscopic two-body scattering at low energy. In the microscopic view, we 
observe that the contribution of the attractive channel exceeds that of the repulsive 
channel at p± = 150 MeV. The effect of the mean field interaction diminishes rapidly 
as the emitted pion recedes from the freeze-out hypersurface. 

In Fig. [9l we show the transverse momentum dependence of the effective pion 
mass at the moment, t = 5 fm/c. On the rightmost panel, we find the enhance- 
ment of the effective pion mass on the backside region which is the effect of the 
repulsive mean field interaction. Owing to the repulsive tttt scattering above the p 
meson mass, the mean field interaction can be repulsive. Although one may consider 



*' With respect to the origin set on the center of the nucleus collision, the detector side indicates 
the source region in the direction of the averaged momentum of the pair. The backside is the other 
half of the source region. 



30 



K. Hattori 



that the attractive irn scattering at low energy dominates the contribution to the 
convolutional integral in Eq. (|4 5j) . the repulsion in the higher regime exceeds the 
attraction owing to the effect of the collective expansion. Because the pion emitted 
in the backside region penetrates into the interior of the source against the collec- 
tive motiorF^I. the center-of-mass energy of the scattering with expanding medium 
is typically large compared with the case emitted on the detector side. Owing to 
the large center-of-mass energy, the repulsive interaction observed on the rightmost 
panel appears at the momentum, p± = 500 MeV, which is less than the mass of p 
meson resonance. 

As mentioned below Eq. (j4- 12|) . we have a negative contribution of the distribu- 
tion function f(t, x,p') owing to the longstanding issue of the Cooper-Frye formula. 
Nevertheless, the enhancement of the effective mass is not an associate artifact of 
this issue. We have confirmed that the negative distribution results in less than 2% 
variation of the self-energy. 

In Fig. [TUJ we show the attenuation factor given by Eq. (|5-10p . which is obtained 
by integrating the imaginary part of the self-energy along the world line. We have 
computed the world line using the classical equation of motion (|5-7p with the real 
part of the self-energy. A set of boundary conditions is imposed by the emission point 
(x, y) and the asymptotic transverse momentum k± indicated below each panel. The 
attenuation factor is displayed using the contour plot with the darker color for the 
stronger attenuation. Each position on the panel (x, y) is the emission point located 
on the freeze-out hypersurface viewed on the transverse plane: therefore, the emission 
at each time is plotted on the same panel. Note that the emission point (x, y) is 
not the apparent spatial coordinate (xf,y') defined by Eq. (|2-64p . but the original 
coordinate before being transformed. 

We find that the attenuation is strong for the pions emitted in the backside 
region, because the path length to escape the interaction range is large for them 
compared with the case emitted in the surface region on the detector side. This 
feature is common to the cases at different momenta. Comparison of the left two 
panels at k± = 100 MeV and k± = 200 MeV shows that the attenuation becomes 
stronger at k± = 200 MeV as the transverse momentum becomes larger. This 
is because the imaginary part of the scattering amplitude shown in Fig. [7] is an 
increasing function of y/s below the p meson mass. On the other hand, we find that 
the profiles of the attenuation factor at k± = 200 MeV and k± = 400 MeV are 
similar to each other almost independent of the transverse momentum. Since the 
imaginary part of the scattering amplitude is dominated by the p meson peak, the 
convolutional integral in Eq. (|4- 5|) converges once the center-of-mass energy y/s goes 
beyond the p meson peak. 

4.2. Modification of the Gaussian radius parameters 

We investigate how and to what extent the dynamical mean field interaction can 
be effective to the observables measured in the ultrarelativistic heavy ion collision. 
Using the phenomenological mean field potential, or the self-energy, obtained in £ 14.14 



Note that we show the effective pion mass for <j> = in Figs. [5] and [5] 
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Fig. 11. Profiles of the distorted images viewed on the apparent coordinate system at k± = 200 
MeV: These panels show the contour plots of the effective distribution function in Eq. (|2-63p . 
The location on each panel is the apparent emission point (xt,y') transformed by Eq. (|5-10[) 
from the original coordinate (x,y). The color becomes lighter as the distribution increases. The 
leftmost panel shows the free case, in which the image is elongated in the outward direction 
owing to the temporal structure of the source, although we have the cylindrical geometry for 
the spatial structure. We show the effects of the real part of the potential on the central panel, 
and the effects of both the real and imaginary parts on the rightmost panel. 
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Fig. 12. Profiles of the distorted images viewed on the apparent coordinate system at k± = 400 
MeV: See the description in Fig. 1111 

we show the distorted profiles of the HBT images, as in the static case examined in 
Sj3j To make these findings more quantitative, we show the effects on the conven- 
tionally analyzed Gaussian radius parameters. The phase shift is evaluated here by 
calculating the relativistic one-particle classical action as described in Appendix A, 
while we have performed nonrelativistic analysis in 

In Figs. [11] and \12\ we show the profiles of the distorted images viewed on 
the apparent coordinate system. Each panel shows the contour plot of the effective 
distribution function defined in Eq. (|2-63p . as in Figs. OandO The results shown in 
Sj3] for the static source is improved here by incorporating the time dependence of the 
freeze-out hypersurface and using the phenomenological mean field interaction. For 
the dynamical source, we find the apparent deformation of the image even in the free 
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Fig. 13. Gaussian radius parameters, R s idc, -Rout and their ratio R out /R B idc ■ The dotted curve 
in each panel indicates the case without the mean field interaction. While the dashed curve 
indicates the effect of only the real part of the mean field interaction, the solid curve indicates the 
effect of both the real and imaginary parts. The effects of the mean field interaction are effective 
below k± < 0.5 GeV. The experimental data measured in Au+Au collision at -^/sjvjv = 200 GeV 
are indicated by the circles and triangles.^ 

cases as displayed on the leftmost panels in each figure. The images are elongated 
in the outward direction owing to the temporal structure of the source, although we 
have the cylindrical geometry for the spatial structure. Only the reflection symmetry, 
with respect to the outward axis, is preserved owing to the preferred orientation in 
which the pion pair is detected. If we analyze a pair carrying averaged momentum 
oriented in the other direction, we find source images obtained by rotating those in 
Figs. H] and [H 

In the free case, the image is more strongly elongated at the higher momentum, 
k± = 400 MeV, since the apparent shift of the spatial coordinate in Eq. (|5T0p is 
proportional to the particle velocity: actually, it is given by the distance, Vk(to — t^), 
where to an d tn are the spatial coordinate dependent emission time and an arbitrarily 
chosen origin of the temporal coordinate, respectively. As we find from the intervals 
of curves, the slope of the image at /cj_ = 400 MeV is much steeper than the slope at 
k± = 200 MeV. It indicates that, at the higher momentum, emissions concentrate on 
the front surface of the source region. This is the effect of the radial flow embedded 
in the thermal distribution on the freeze-out hypersurface, which has been observed 



Distortion of HBT Images 



33 



in the hydrodynamical picture P^* Since the magnitude of the radial flow is large in 
the surface region, energetic pions are likely to be emitted in the thin surface region 
on the detector side, owing to the strong boost of their momenta. 

The remaining four panels show the effects of the mean field interaction. The 
central panel in each figure shows the effect of the real part, and the rightmost 
panel shows the effects of both the real and imaginary parts. The effect of the real 
part stretches the image in the sideward direction at k± = 200 MeV, and is less 
effective at k± = 400 MeV. Comparing with what we observed in Figs. [2] and [3] 
with a schematic potential, we find that these effects are caused by the attractive 
mean field interaction, which is obtained in §4.11 Comparison of the rightmost 
panel with the central panel shows that the slope along the outward axis is steeper 
on the rightmost panel, which indicates a large distribution on the detector side. 
The effect of the absorption cuts off the emission in the backside region, and the 
surviving particle tends to come from the region on the detector side. As we have 
mentioned below Eq. (|2-63p . the normalization of the effective distribution function 
is preserved to be one at each momentum even in the presence of the absorption, 
owing to the same absorptive effect on the single-particle spectrum. Therefore, the 
effective distribution function shows the relative distribution of the emission points 
of the surviving particles. 

To investigate these effects more quantitatively, we show the Gaussian parame- 
ters i?out) R S ide an d their ratio -R ut/-R s idc i n Fig- CG3 These analyses are based on 
the formulae given in Eqs. (|2-33j) . (|2-57p and (|2-58p . The horizontal axis in each plot 
is the magnitude of the asymptotic transverse momentum k± carried by a pion pair. 
The dotted curves display the free cases. They are the Gaussian fitting parameters 
of the distribution at the "freeze-out" specified by Eqs. (|4-8p and (|4-9p . and qualita- 
tively reproduce the results obtained by the hydrodynamical simulations.^ While 
the Gaussian parameter for the sideward extension, i? s id e > measures the bare spatial 
extension, that of the outward extension, R ou t, is enhanced by the temporal struc- 
ture. They have exactly the same value at k±_ = 0, since we have no coordinate shift 
in Eq. (|2-33|) for vanishing velocity, v± = 0. The filled triangles and circles indicate 
the data measured in the central Au+Au collision at y / ijv r /v = 200 GeVpD 1 '^ 1 

We show the effects of the mean field interaction with dashed curves and solid 
ones. The dashed curves indicate the effect of the real part of the mean field poten- 
tial, and the solid ones the cooperative effects of the real and imaginary parts. We 
find that the sideward radius i? s ide is increased by the effect of the real part, and that 
it acts more effectively in the low-momentum regime. This behavior qualitatively 
agrees with what we observed in Figs. [11] and PT2l and also in Figs. [2] and [3] with a 
schematic attractive potential. We find that the deviation from the data in the side- 
ward radius is improved by the phenomenological mean field interaction. However, 
the pion density in the vicinity of the freeze-out hypersurface is not sufficiently large 
to solve the puzzle completely. 

The absorption unexpectedly acts to reduce the sideward radius -R s ide> contrary 
to the findings in Figs. [2] and [3J There, we have argued that the extension in the 
sideward is stretched by cutting off the emissions in the interior and backside region 
of the source. This conflict is due to two main reasons. One reason is that our mean 
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field model incorporates the interactions only among the emitted pions, without 
the mean field interaction by the thermalized matter. If we incorporate the latter 
interaction, the absorption by the matter acts to increase the sideward radius as we 
have already discussed. This effect could be examined in the consistent description 
of the mean field interaction throughout the dynamics of the matter. 

The other reason is found in the profiles of the attenuation factor shown in Fig. 
[TOl Focusing on the variation of the magnitude along the line parallel to the outward 
axis, we find that the attenuation factor becomes larger as the distance from the 
sideward axis increases. This behavior is clearly seen at k± = 100 MeV, and is less 
pronounced at k± = 200 MeV and k± = 400 MeV. Such a profile of the attenuation 
factor leaves a large distribution in the interior of the source compared with that 
in the peripheral region, and it obviously results in the reduction of the sideward 
radius. Recall that the imaginary part of the convolutional integral in Eq. (j4-5|) is 
dominantly contributed by the p meson peak, and that it has little contribution from 
the high-energy scatterings beyond the range of the peak. Therefore, the integral 
monotonically increases as more energetic scatterings contribute to it, unless the 
scattering energy goes beyond the range of the peak. Because the energetic particle 
is likely to be emitted in the surface region owing to the effect of the radial flow, 
the integral has much contribution from the energetic scatterings in that region. 
These contributions cause the observed large attenuation factor in the peripheral of 
the source region at k± = 100 MeV. As the momentum of the escaping pion, k±, 
increases, the attenuation factor becomes large in the interior of the source, while 
it is convergent in the surface region with the scattering energy going beyond the p 
meson peak. Thus, the attenuation factor depends on the sideward coordinate much 
weakly at at k± = 200 MeV and k± = 400 MeV. Note that the attenuation factor 
of the pion emitted in the front surface region on the detector side is not large, since 
it drifts for the collective motion of the medium, and the center-of-mass energy is 
typically small. 

The effect of the absorption slightly reduces the outward radius i? ut- This 
behavior qualitatively agrees with our intuition and the clear observations in Figs. [2] 
and [3j that is, the emissions in the backside region are cut off by the absorption. The 
magnitude of the effect is, however, not so strong. It is largely because the mean 
field interaction by the matter is not incorporated, as mentioned above. Further 
investigations are needed to examine the opacity due to the mean field interaction. 
Owing to the interaction with the matter, the effect of the absorption would be 
enhanced and act to reduce the deviation in the outward radius. 

Owing to the attractive interaction, the ratio, -R ut/-R s ide> is improved by 15% 
below 300 MeV. The magnitude of the effect is, however, quantitatively insufficient 
to resolve the deviation completely. By the hydrodynamical simulations, the ratio is 
overestimated owing to the overestimation of i? ou t , and the underestimation of R sl dc , 
in contrast. Thus, the effect of improving the deviation should act on the radii in a 
different way. 

For the sake of completeness, we show the modification of the transverse pion 
spectrum at the central rapidity. In Fig. [21 the dotted line displays the spectrum in 
the free case, which freezes abruptly at the hypersurface. The circles and triangles 
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Fig. 14. Modification of the transverse pion spectrum at the central rapidity: The dotted line 
indicates the free case, which is fixed on the freeze-out hypersurface, and the thin solid line the 
effects of the real part of the mean field potential. The thick solid line shows the effects of both 
the real and imaginary parts of the mean field potential. We show these effects in the linear and 
logarithmic scales, respectively. The circles and triangles exhibit the spectrum of the positive 
pion measured in Au+Au central collision at ^/snn = 200 (GeV) in the RHIC experiment 
The negative pion is measured in almost the same amount. 



indicate the spectrum of the positive pion measured in Au+Au central collision at 
V^v" = 200 GeV in the RHIC experiment P>® The negative pion is measured in 
almost the same amount. In the linear scale, we find that the effect of the attractive 
interaction acts to increase the distribution in the low-momentum regime, and to 
decrease that in the higher regime around 0.3 to 0.4 GeV. These effects show a 
deceleration of the emitted pions by the attractive mean field interaction. However, 
the effect of the real part is so slight that we find little deviation in the logarithmic 
scale. The distribution in the regime higher than 0.5 GeV is barely affected by the 
real part of the mean field potential, because the shift of the momentum is much 
smaller than the magnitude of the momentum itself. The thick solid line indicates 
the effects of both the real and imaginary parts of the mean field potential. The 
absorptive effect reduces the amount of pions below 0.6 GeV, which is mainly due 
to the formation of p meson resonance, tttt — > p. Because the inverse reaction is not 
incorporated in our model as we have mentioned below Eq. (|2-63p . the microscopic 
detailed balance is not maintained 1*^1 To refine the model, the finite lifetime of the 
resonance should be taken into account. Considering the decay process, p — > tttt, 
which has the almost 100% branching ratio, the pion number is conserved within the 
binary tttt collision. 

Finally, we comment on the validity of the semiclassical method used in this 
analysis. If the classical action S c satisfies a condition 



d 2 S c 




dS c 


dx 2 


< 


dx 



(4-17) 



*' Nonetheless, the normalization of the effective distribution function in Eq. (|2-63p is preserved 
even in the presence of the absorptive effect, as mentioned there. 
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we have Hamilton-Jacobi equation as the classical limit of the Schrodinger equation, 
and the approximation works well. The condition ()4- 17|) is rewritten as 



A 



y/2m(E - V(x)) 



< 2(E-V(x)) 



dV{x) 



dx 



(4-18) 



which indicates that the spatial profile of the potential V(x) should not vary rapidly 
over a wavelength A, and that the condition tends to be satisfied in the high-energy 
limit. Our computation is based on the semiclassical approximation applied to the 
four-dimensional pseudo-Schrodinger equation (|5-2|) . in which "energy", "potential 
" and "mass" are given by E = —m?/2, V{x ll ,p) = 77(x M ,p)/2 and fh = 1, respec- 
tively. Thus, we have a condition to justify the approximation as an extension of 
Eq. (|4T8p . Using the physical quantities, the condition is represented by 



m cS \ -3 1 
m J 



2m? 



dx^ 



< 1 



(4-19) 



where m e g is the effective mass defined in Eq. ()4T6p . Evaluating the left-hand side 
of Eq. (|4-19|) . we obtain a sufficiently small value not larger than 0.05. Therefore, 
the semiclassical method employed in our analysis works well to evaluate the phase 
shift of the amplitude caused by the mean field interaction. 



§5. Concluding remarks 



In this work, we have examined the effects of a final state interaction on HBT 
interferometry in the ultrarelativistic heavy ion collision. This study is originally 
motivated by the so-called "RHIC HBT puzzle", that is, an observation that the 
theoretical estimates of the HBT Gaussian radii show a systematic deviation from 
the experimental values measured in the RHIC experiments. To find a solution to this 
puzzle that has been addressed so far, we have studied how the mean field interaction 
ignored in the conventional hydrodynamical modeling may cause the distortion of 
the source images. 

First, in ^21 we have examined how the effects of the mean field interaction are 
incorporated in the framework of HBT interferometry, beginning with the definition 
of the correlation function in terms of the density matrix. We found that the differ- 
ence in the phase shifts imposed on each one-body amplitude of a pair induces the 
coordinate transformation in (|2-59p . It maps the original space-time coordinate onto 
the apparent one. The original grid, mapped onto the apparent coordinate system, 
is deformed owing to the nonlinear nature of the transform, as shown in Fig. [TJ 
Then, the images viewed on the apparent coordinate system exhibit the distortions. 
We have shown that this is the geometrical interpretation of the distortion of the 
HBT images.^ 

In ^31 we have shown that the images are distorted contrastingly depending on 
whether the interaction is repulsive or attractive: a repulsion acts to shrink the 
sideward extension and an attraction acts to stretch it, in contrast. The effect of the 
absorption cuts off the emissions in the interior and backside region of the source. 
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This effect stretches the images in the sideward direction, which is somewhat similar 
to the effect of attractive interaction. 

In E j4.ll a phenomenological model of the mean field interaction is constructed 
on the basis of the elastic forward tttt scattering amplitude and the pion distribution 
function after the "freeze-out". Owing to the dominant contribution of the p-wave 
scattering, we found the attractive mean field interaction in the low-momentum 
regime in some hundreds MeV. In this regime, the mean field interaction acts on 
the images effectively. In the high-momentum regime, the attraction turns into 
repulsion due to the repulsive regime in the p-wave scattering beyond the p meson 
mass. However, we did not find a considerable effect of the repulsive interaction. 

Using the phenomenological self-energy obtained in £ )4.1l we have computed the 
modification of the amplitude of an emitted pion, that is, the phase shift and atten- 
uation factor. We have employed the semiclassical method given in Appendix A. We 
examined how and to what extent the modification results in the Gaussian radius 
parameters. We found that the effect of the obtained attractive interaction actu- 
ally increases the sideward radius -R s ide as implied in the preceding works EliQ3J2D 
and that it improves the deviations in i? S ide and the ratio, R out /R s idc ■ However, 
the magnitude of the effect is not sufficiently strong to resolve the deviations com- 
pletely, since the pion density is not sufficiently large in the exterior of the freeze-out 
hypersurface owing to the evaporation. 

We briefly comment on the related works. It has recently been suggested that 
some upgrades of the hydrodynamical model lead to considerable improvements of 
the deviations owing to the cooperative effects of the upgrades.^ 1 We note that the 
modification of R s ide obtained in this work is on the same order of their individual 
effects. The mean field interaction in the matter, not included in our model, would 
act to vary the apparent radii more efficiently, since the density is larger than that in 
the very last stage. This possibility should be examined with the consistent descrip- 
tion of the mean field interaction throughout the hadron phased Alternatively, if 
we attempt to treat the effects of the mean field interaction in the hydrodynamic 
picture, we could study it with the equation of state. Its effects possibly reflect in 
the profile of the freeze-out hypersurface. 
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Appendix A 

We examine the modification of the amplitude incorporating the relativity. The 
amplitude of the pseudo scalar particle obeys the Klein-Gordon equation 

(n + m 2 + n(x))if k (x) = , (5-1) 

where the coordinate x denotes the four- vector x^ = (t,x), and the box stands 
for the derivative operator, □ = d 2 — d 2 . The effect of the mean field interaction 
is incorporated in the self-energy which is the mass shift of the pion mass, 

m, in the vacuum. The modification of the dispersion relation in the pion cloud 
is phenomenologically studied in §4.1 1 In this appendix, we examine a method of 
obtaining the semiclassical form of the amplitude. 

A parameter called the proper time enables us to adopt the semiclassical approx- 
imation to the Klein-Gordon equation (|5T|) employing the standard prescription in 
quantum mechanics] E2DJ37} Tj s j n g p rQ per time A , the Klein-Gordon equation 
(|5Tp reads 

° + 5j!l)sp fc (x,A) = i^ fc (x ) A) ) 

m 2 ( 5 ' 2 J 

!Pfc(s,A) = e % T x <p k (x) , 



where <Pfc(a;, A) is the amplitude that obeys the pseudo-Schrodinger equation (|5-2p 
with the unit mass, and the pseudo potential is given by the self-energy, 77/2. The pa- 
rameter A plays a role of "time" in nonrelativistic quantum mechanics, and !Pfe(x, A) 
is the eigenfunction of "energy" given by —m 2 /2. We deduce that the pseudo- 
Schrodinger equation (|5-2p is obtained from the Hamiltonian in the four dimensions, 

H(*,P) = -J + E P-, (5-3) 

where we take a convention of the Minkowski metric, p 2 = Pq — p 2 . 

Owing to the argument in Q2.2\ the transition amplitude ip(x' , t';x,t) from (t, x) 
to (t',x') provides the amplitude tpk(%) of a particle emitted at (t, x) to detect with 
the asymptotic momentum k, as long as the classical trajectory is uniquely specified 
by x and k. Assuming that the transition amplitude ip(x' ,t';x,t) obeys the Klein- 
Gordon equation (|5-ip . the technique with the proper time enables us to write 

<p{x', t'; x, t) = e- in ^ x (x'\e- iRx \x) . (54) 

On the right-hand side, the transition amplitude, from x^ to x'^ in the "time" interval 
A, obeys the pseudo-Schrodinger equation (|5-2p . Using the conventional prescription, 
we obtain the path integral form of the amplitude: 

(x'\e- iSx \x) = !®x fape'ti** {-p^-h^p)) , (5.5) 



where the acute on the coordinate denotes the derivative with respect to the proper 
time as 

dt dx 



d\ ' r d\ 



-p^x^ = -p — +p 
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We evaluate the right-hand side of Eq. (|5-5|) adopting the stationary phase 
approximation. The variation of the real part, 5(p )l x^ + ReH(x,p)), leads to the 
canonical equations of motion 

C £/i _ dS r 

where H r is the real part of the Hamiltonian, H r {x,p) = Re H(x,p). The proper time 
parameterizes the world line of the classical motion. Using the temporal component 
of Eq. (|5-6p . we eliminate the proper time in the spatial components and parametrize 
the world line using the physical time, t. Substituting the Hamiltonian (|5-3p . the 
spatial components of Eq. (|5-6p are expressed as 



x 



p 

E,, 



v 



■ n _ i an r ( x ,t) O 7 ) 

[ P ~ 2E P dx ' 

where II r (x) is the real part of the self-energy, II r (x) = Rei7(x), and the dot 
indicates the derivative with respect to time, t. The zeroth component of p^ is 
given by E p = \/p 2 + m 2 + 7T r . The set of canonical equations (|5-7p determines the 
classical trajectory under the existence of the mean field interaction. 

During the classical motion, the "energy" given by —m 2 /2 is conserved, be- 
cause the Hamiltonian (|5-3p does not depend on the proper time. Then, we have a 
dispersion relation, 

- p 2 n r (x) m 2 

H r (x,p) = = , 

rv if) 2 2 2 

where the Hamiltonian is regarded as c-number under the semiclassical approxima- 
tion. Substituting the above expression and Eq. (|5-5p for Eq. (|5-4p . we obtain the 
transition amplitude 

<p{x', t'; x, t) = ^ e ifp d*-jE P dt) ^ ( 5 .g) 



where the integrals are performed along the world line determined with Eq. ((5Z 
Owing to the change in the integral variable using Eq. (|5-7p . we find a simple 
extension to the relativistic form, 

<p(x',t';x,t) = ^e iS{x '' t '- x > t) , 

rt' 

S(x', t'; x,t) = - J mcffVl - x 2 dt' . (5-9) 

The effect of the mean field interaction is incorporated in the effective mass, m e g = 
appearing in the relativistic one-particle action S(x' ,t';x,t). Taking 
the imaginary part of the Hamiltonian in Eq. (|5-5p . the prefactor stf associated 
with the classical trajectory is provided by integrating the imaginary part of the 
self-energy along the world line, 

srf = exp I — -Im [n(x)] dt , (5-10) 
J 2£U 
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where we have used the temporal component of Eq. (|5 6j) to take time, t, as the 
integral variable. 

As mentioned above, the semiclassical evaluation of <p(x', x) provides the desired 
amplitude cpk(x). As the four-dimensional analogue of Eq. (|2-38[) . we have 

ip k (x) = Jdx' e ikx ' ip(a/,x) , (5-11) 

where x is the coordinate four vector, and k is the asymptotic momentum that 
satisfies the on-shell condition, E 2 , = k 2 + m?. If we choose sufficiently a large time 
t\ and spatial coordinate x\ far outside the range of the mean field potential, the 
relativistic classical action in Eq. (|5-9[) is decomposed into two parts as in Eq. (|2-39p . 

S(x', x) = S(x', x\) + S(x±,x) , 
= k{x — x\) + S{x\,x) , 

where S(x',x\) is the classical action of the free motion from x\ to x' . Inserting the 
above relation into Eq. (|5-lip and performing the integral, the amplitude is given 
by 

<p k (x) = st e ~ikx + iSS(x u x) ( 5 . 12 ) 

5S(xi,x) = S(xi,x) — k(x\ — x) , (5 - 13) 

where 8S(x\,x) is the phase shift caused in the range of the mean field potential. 
Using the formulae in Eqs. ()5-12p and (j5- 13fl . we obtain the amplitude from the 
computation of the relativistic one-particle action, the integral of which is performed 
along the world line determined using the canonical equations (|5-7p . 

We remark that the dynamical mean field potential, or phenomenological self- 
energy, studied in ^4.11 depends on the momentum as well as the space-time coor- 
dinate. Then, the canonical momentum appearing in the Hamiltonian should be 
distinguished from the kinetic momentum defined by the velocity. However, we ne- 
glect their difference which is given by the derivative of the self-energy with respect 
to the momentum, since it is negligibly small in our model compared with the mag- 
nitude of the momentum in some hundreds of MeV. 
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